GIbPSs: a toolkit for fast and accurate analyses of genotyping-by-sequencing data without a reference genome.

Genotyping-by-sequencing (GBS) and related methods are increasingly used for studies of non-model organisms from population genetic to phylogenetic scales. We present GIbPSs, a new genotyping toolkit for the analysis of data from various protocols such as RAD, double-digest RAD, GBS, and two-enzyme GBS without a reference genome. GIbPSs can handle paired-end GBS data and is able to assign reads from both strands of a restriction fragment to the same locus. GIbPSs is most suitable for population genetic and phylogeographic analyses. It avoids genotyping errors due to indel variation by identifying and discarding affected loci. GIbPSs creates a genotype database that offers rich functionality for data filtering and export in numerous formats. We performed comparative analyses of simulated and real GBS data with GIbPSs and another program, pyRAD. This program accounts for indel variation by aligning homologous sequences. GIbPSs performed better than pyRAD in several aspects. It required much less computation time and displayed higher genotyping accuracy. GIbPSs retained smaller numbers of loci overall in analyses of real GBS data. It nevertheless delivered more complete genotype matrices with greater locus overlap between individuals and greater numbers of loci sampled in all individuals.

Hybridization between mouse lemurs in an ecological transition zone in southern Madagascar.

Hybrid zones in ecotones can be useful model systems for the study of evolutionary processes that shape the distribution and discreteness of species. Such studies could be important for an improved understanding of the complex biogeography of Madagascar, which is renowned for its outstanding degree of small-scale endemism. Certain forest remnants in central Madagascar indicate that transitional corridors across the island could have connected microendemics in different forest types in the past. Evolutionary processes in such corridors are difficult to study because most of these corridors have disappeared due to deforestation in central Madagascar. We studied a hybrid zone in one of the few remaining ecotonal corridors between dry and humid forests in Madagascar, which connects two species of mouse lemurs, Microcebus griseorufus in dry spiny forest and Microcebus murinus in humid littoral forest. We sampled 162 mouse lemurs at nine sites across this boundary. Morphometric analyses revealed intermediate morphotypes of many individuals in transitional habitat. Bayesian clustering of microsatellite genotypes and assignment tests yielded evidence for a mixed ancestry of mouse lemurs in the ecotone, where we also observed significant linkage disequilibria and heterozygote deficiency. In contrast to these observations, mitochondrial haplotypes displayed a sharply delimited boundary at the eastern edge of spiny forest, which was noncoincident with the signals from microsatellite data. Among several alternative scenarios, we propose asymmetric nuclear introgression due to male-biased dispersal, divergent environmental selection, and an expansion of dry spiny forest in the course of aridification as a probable explanation of our observations.